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Abstract 

Statistical properties of eigenvectors in non-Hermitian random matrix ensem- 
bles are discussed, with an emphasis on correlations between left and right 
eigenvectors. Two approaches are described. One is an exact calculation for 
Ginibre's ensemble, in which each matrix element is an independent, identi- 
cally distributed Gaussian complex random variable. The other is a simpler 
calculation using N~'^ as an expansion parameter, where A'^ is the rank of 

'^ . the random matrix: this is applied to Girko's ensemble. Consequences of 

^ , eigenvector correlations which may be of physical importance in applications 

O \ are also discussed. It is shown that eigenvalues are much more sensitive to 

perturbations than in the corresponding Hermitian random matrix ensem- 
bles. It is also shown that, in problems with time-evolution governed by a 

Q>^ . non-Hermitian random matrix, transients are controlled by eigenvector cor- 

C^^ , relations. 
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I. INTRODUCTION 

Hermitian random matrices have been very successfully used to model Hamiltonian opera- 
tors of closed quantum systems |jl| . In many cases, this has lead to a quantitative description 
of features such as spectral fluctuations in classically chaotic quantum systems and in disor- 
dered quantum systems in the metallic regime 0. Within this approach it is also possible 
to describe the statistical properties of wave functions and matrix elements in such systems. 
Random matrices are also of great importance in many other areas of physics in which they 
are not constrained to be Hermitian 0,3- These include: the dynamics of neural networks 
^, the quantum mechanics of open systems [^, classical diffusion in random media ^ 
and in population biology p, and modelling the statistical properties of flux lines in su- 
perconductors with columnar disorder P|-p^. Recently, in connection with these problems, 
spectral properties of non-Hermitian random matrices and operators have been studied in 
great detail (see for instance [§|,[7|,|,|T3-|I|]). 
In the context of fluid dynamics it is well known ||T9|^T| that systems governed by non- 



Hermitian evolution operators exhibit striking features. First, such systems are particularly 
sensitive to perturbations. Second, these systems can exhibit pseudo-resonances at which the 
system reacts strongly to an external perturbation although the excitation frequency is not 
close to any of the frequencies of the internal modes. Third, non-Hermiticity can give rise to 
interesting transient features in time evolution. Such features cannot be understood solely 
in terms of the spectrum of the evolution operator. While the eigenvalues of the evolution 
operator determine the long-time behavior, transients and sensitivity to perturbations, in 
particular, are determined by the properties of the corresponding eigenvectors. 
In this paper we quantify the statistical properties of the eigenvectors of random non- 
Hermitian matrices and examine to what extent enhanced sensitivity to perturbations and 
transients in time evolution are present in random systems described by non-Hermitian op- 
erators, such as Fokker-Planck operators or projected Hamiltonians [||. We report on 
results for two ensembles of random matrices, namely Ginibre's ensemble and Girko's 
ensemble 0]. There are several reasons for studying these ensembles. First, characterizing 
the statistical properties of eigenvectors in these cases is by itself a problem of consider- 
able interest: we show that left and right eigenvectors exhibit striking correlations, which 
depend strongly on where in the spectrum the corresponding eigenvalues lie. Second, non- 
Hermitian operators (such as the Fokker-Planck operator governing classical diffusion in 
a random velocity field 0) may be represented, in a finite system and in an appropriate 
basis, as random matrices. In general, their matrix elements exhibit certain structures and 
are much less uniform than the matrices from the ensembles we investigate. Nevertheless, 
experience of universality in random Hermitian problems gives reason to hope that random 
matrix ensembles will provide a widely-applicable guide to behavior. Third, the high sym- 
metry of Ginibre's ensemble (the matrix elements are independently Gaussian distributed) 
allows for an exact calculation which we present in detail. Separately, we develop an alterna- 
tive, more general and simpler approach to calculations, based on a perturbative evaluation 
of ensemble- averaged resolvents, using A^~^ as the expansion parameter, where N is the rank 
of the matrix. We apply this to Girko's ensemble, and also assess its validity in Ginibre's 
ensemble by comparison with exact results. Such approximate methods are particularly 
important because they are easily extended to the more general ensembles discussed in MM. 



The remainder of this article is organized as follows. In section H we discuss the formulation 
of the problem: we define the ensembles of random non-Hermitian matrices studied in sub- 
sequent sections, define the densities of eigenvector overlaps that will be the main objects 
of study in this paper, and the corresponding Green functions. This section establishes the 
notation used subsequently. In section |T| we show how to derive exact results for the sta- 
tistical properties of eigenvectors in Ginibre's ensemble of non-Hermitian random matrices 
of arbitrary matrix dimensions. We also discuss simplifications which arise in the large N 



limit. In section |I^ we summarize our approximate calculations, applying them to Girko's 
ensemble. The eigenvector correlators are calculated in terms of two-point functions which 
are obtained within the self-consistent Born approximation. This approach is appropriate 
in the limit of large N, and under certain additional assumptions which are discussed in 
this section. As a special case, we check results found with this method for Ginibre's en- 



semble against the exact results of section |T|. The results obtained in sections |T| and |V 



are summarized and discussed in section |V|. In section [V^ we examine some consequences 
of eigenvector correlations that are likely to be important in physical applications. These 
are: the extreme sensitivity of eigenvalues to perturbations; time evolution governed by 
non-Hermitian random matrices; and the nature of correlations between individual eigen- 
vector components in non-Hermitian random matrix ensembles. Finally, we summarize our 
conclusions in section |VI1| . An outline of some of these results has been published previously 



in a shorter communication 22 



II. FORMULATION OF THE PROBLEM 

A. Ensembles of non-Hermitian matrices 

In the recent literature, a large number of different ensembles of non-Hermitian random 
matrices and operators have been discussed [^0,0,||. In the following, we restrict ourselves 
to Ginibre's and Girko's ensembles [Q of non-Hermitian random matrices. 
Ginibre introduced an ensemble of random N x N matrices J which have complex elements 
Jki with independently distributed real and imaginary parts J^^ and J^';. The ensemble is 
defined by the measure [^ (see also |Q) 

P{J) d J (X exp ( - 1 Tr [Jjt]) J] dJ^.d J^;, . (1) 

kl 

Thus (Jki) = and {JhJm) = o"^- Here (■ ■ ■) denote ensemble averages and the overbar 
indicates complex conjugation. The parameter a"^ controls the density of eigenvalues: in the 
limit N — > oo, the ensemble- averaged density (per unit area) is l/vro"^ within a disk in the 
complex plane, centered on the origin and of radius yNa. Two different conventions are in 
use for the value of a^. The choice a^ = 1 (as for instance in |jl|,|^) results in a fixed density 
as A^ -^ CX3. Alternatively, the choice a^ = A^~^ results in a fixed support for the eigenvalue 
density as A^ — >• cx). 
Girko has considered the following generalization of Ginibre's ensemble, 

P(J)dJcxexp('-1^47lTr[JJ^-rReJJ]') \{dJlid4,, (2) 

\ ^ '^ J kl 



with — 1 < r < 1. In this ensemble the non-zero cumulants are 

{Jkl Jkl) = cr^ , {Jkl Jlk) = TO-^ • (3) 

For r = 0, Ginibre's ensemble is recovered; the case r = 1 corresponds to Dyson's Gaus- 
sian Unitary Ensemble |jl|, while r = — 1 describes an ensemble of complex anti-Hermitian 
matrices. 

B. Densities of left and right eigenvectors 

The eigenvalues, Aq, and left and right eigenvectors, (Lq,| and \Ra), of the matrix J satisfy 

J\Ra) = Aq, \Ra) , (4) 

In general, the eigenvalues are complex numbers Aq, = A^ + ^^a- Except for a set of measure 
zero, they are non-degenerate. In this case the eigenvectors form two complete, biorthogonal 
basis sets with the normalization 

{La\Rp) = 5ap. (5) 

The closure relation is 

Y,\La){Ra\ = l. (6) 

a 

We denote the Hermitian conjugates of (Lq,| and \Rp) by |Lq,) and (-R/3I, so that, for example, 
\La) satisfies J^\La) = A*|Lq,). Left and right eigenvectors are generally not orthogonal 
amongst themselves. On the contrary, scalar products can vary significantly. This can 
have important physical implications. For instance, it is well-known that non-orthogonality 
of eigenvectors can have an important bearing of time evolution in systems governed by 



non-normal operators |21 



In the following we consider statistical properties of scalar products of eigenvectors in ensem- 
bles of random non-normal operators. We note that Eqs. (^ and (|]) allow for the following 
scale transformation 

\Ra} — ^ Ca \Ra} , \^ci\ — ^ \^a\ Cq, (7) 

with arbitrary complex numbers Ca'- we study only such combinations of eigenvectors as are 
invariant under this scale transformation. The simplest such combination of two eigenvectors 
is trivial [see Eq. (^]. We hence consider the combination 

OaP = {La\Lp) {Rp\Ra) . (8) 

We calculate the mean value and discuss the distribution function of this overlap matrix. 
Note that completeness implies the sum rule 

E^"/3 = i- (9) 



It is convenient to define local averages of diagonal and off-diagonal elements of Oap-, 

0{z) = la''Y,0aa5{z - XoS) , (10) 



0(zi, Z2) = {a^ Y. ^-/3 Kzi - K) Kz2 - A^)) . (11) 

Here, z = x + iy is a complex number with real and imaginary parts x and y and 5{z) 
denotes a delta-function in both coordinates. Correspondingly, the density of states and the 
two point function are defined as 

d{z) = (a^Y.^{z - X^)) (12) 

a 

R2{zi, Z2) = (a^ Y, K^i - K) S{Z2 -Xp)). (13) 

In order to characterize the overlap matrix using Green functions, it is convenient to intro- 
duce the density 

D{zi, Z2) = (a^ Yl ^-/3 '^('^1 - ^") ^(-^2 - A/3) ) , (14) 

Q,/3 

which can be expressed in terms of 0(2:1) and 0(2:1, Z2) as 

D{zu Z2) = 0(^1) 5(^1 - Z2) + 0(zi, Z2) . (15) 

Thus, information on the diagonal overlap matrix elements may be extracted from the 
singular part of D{zi,Z2)- The smooth part conveys information on the off-diagonal overlap 
matrix elements. 
Finally we note that the sum rule (|^) implies the constraint for the density -D(2;i, Z2) 

fd^Z2D{z,,Z2)=d{z,) (16) 

where d{zi) is the density of states [Eq. (|I^)]. 

C. Green functions and spectral densities 

We shall make use of the fact that the densities d{z) and D{zi, Z2) may be expressed in terms 
of ensemble averages of resolvents {z — J)~^ and products of resolvents (2:1 — J)~^ (^2 ~ J^)~^. 
The density of states d{z), for example, by means of the relation 

S(z) = i| i (17) 

IT OZ Z 



may be expressed as ||7|JT5 



d{z) = ^-^{TT[{z-jr']). (18) 



In Eqs. (0) and (|T| 



9__lfd_ .d_\ 

Eq. (p!?!) replaces the relation 6{E) = vr^^Ini {E — iO^)^^ which is applicable in problems for 
which the Green functions are analytic in the upper and lower complex half-planes. Here, 
this is not the case. 
Similarly, the density D{zi,Z2) can be obtained from a two-point function 

Di^i^ ^^) = 5 J; ^ <Tr [(., - jr irz2 - J^'] ) . (20) 

This is most easily seen from the spectral representation of the resolvent: [z — J)~^ = 
Sa l-^a) (^ ~ ^a)~^ {La\- Wc show in sectiou |V| how the averages in Eqs. (|T8|) and ( pOD can 
be calculated perturbatively, using an expansion in powers of N~^. 



III. GINIBRE'S ENSEMBLE 

As pointed out in the introduction, Ginibre's ensemble is a special case of Girko's family 
of ensembles of non-Hermitian matrices. It is obtained by setting r = in Eq. (|^) and is 
thus the ensemble of complex matrices with independent, Gaussian distributed elements. In 
this special case we are able to provide an exact calculation of the eigenvector correlators 
introduced in section 0. 

A. Density-of-states and eigenvalue correlations 

Eigenvalue correlations for the ensemble (|l|) were first studied by Ginibre @. The joint 
probability distribution of the eigenvalues is 

P^(Ai, . . . , A^) = C^ n I^M - A.I' exp ( - ^ 5^ |AJ2) (21) 

with normalization C^ = {N\ Y\j=q ^ J' o"^-'"^^)"^. The eigenvalue density and the two-point 
function, d{z) and R2{zi, Z2) [Eqs. (p!^,([T3|)], may be calculated by averaging 6{z — Ai) and 
(A^ — 1) 6{zi — Ai) 6{z2 — A2) with the weight P/y. In the following we demonstrate briefly 
a way of performing the corresponding A-integrals which can be readily generalized to deal 



with the integrals that arise in the calculation of the eigenvector correlations [Sec. |IIIB4 
Making use of the fact that 



P^(Ai, . . . , A^) = (vr AT! a^^)-' exp (- ^) J] J^i - A„|2p^_i(A2, . . . , A^v) (22) 

we have 



m=2 



d{z) = Na^7rN\(r^^)-^ exp 
This can be written as 



I 12 /• ^ 

(-^) d''\2---d''\NPN-liX2,--.,XN)ll\z-\m\''- (23) 

■^ m=2 



d{z) = Na^nNla^' exp 



(j^ 



det 



doQ dQi 
diQ 



dN~3N~2 
dN-2N-3 dN-2N-2 J 



(24) 



with d,j = {nj\a^^+^y^ Jd^XXx' \z - Ap exp(-(T-2 |A|2). Denoting the (A^-1) x (iV-1) 
determinant in Eq. ( plf ) by -Dat-i we derive the recursion relation 

(25) 



Dk+i = {cr-^ \z\^ + k + l)Dk - a-2 \z\HDk-i. 
Using Di = I + cr"^ \z\'^ and D2 = 2 + 2cr~^ |zp + a"'' |2;|'*, we thus obtain 



N-l 






(26) 



which corresponds to Eq. (51.1.32) in [^. In the hmit of large N, with a"^ = N ^, the 
density of states is 



d{z) 



71 ^ for |2;| < 1 , 
otherwise . 



Similarly, we obtain for the two-point function 

R2iz,,Z2) = Na^ (vr^iVlae)-! \z, - z^l' expi-a'^ \z,\^ - a-^\z2\') 



(27) 



(28) 



with 



X det 



/oo /oi /02 

/lO /ll 

/20 



fN-5N-3 

fN-4:N-4 In-AN-S 

fN-3N-5 fN-3N~4 fN-3N-3 



Al 



/,, = (vr (j + 1)! a2^+6)-i / d'xXx' \z, - Ap |^2 - Ap exp ( - ^) 



'u — V" VJ I '-J' ^ ) ^ ^ /\ /\ /\ 1^1 ^v| 1^2 ^'1 exp y —) . v-^'^j 

As before, we derive a recursion relation for the (A^ — 2) x (A^ — 2) determinant in Eq. (^B|). 
This recursion relation simplifies considerably when Z2 = 0. Denoting the determinant in 

Eq. (pSD by Fjv_2, we have with zi = z 



8 



Fk+i = (a-2 \z\^ + k + 2)Fk - a-2 \z\^ {k + 1) F,,_i . (30) 

In this way we obtain, with Fi = 2 + cr^^ |zp and F2 = 6 + 3cr~^ |zp + a~^ \z\^ 

1 



i?2(^,0) 



TT^cr^ 



exp 



a^ 



i-E 



A'-l .,,2, 



«=0 



/!a2' 



(31) 



which (for o"^ = 1) is equivalent to (15.1.30) in [|I| with n = 2 and Z2 = 0. Moreover, in the 
hmit of large A^, with cr^ = N~^, one finds that for zi 7^ Z2 and \zi\, \z2\ < 1, the two-point 
function is constant illl 



N~'R2{ZuZ2) 



71 2 for Zi 7^ Z2 and \zi\, \z2\ < 1, 
otherwise. 



(32) 



B. Eigenvector correlations 

In this section we show how to obtain expressions for correlations of eigenvectors in Ginibre's 
ensemble. We start from ( PUJ ) and (|TTp, perform calculations for general o"^ but set a^ = A^~^ 
in the final results [see Eqs. ( |68D and (|73| ) below]. 



1. Change of basis 
Since the fluctuations of the eigenvectors and those of the eigenvalues are correlated, it is 



convenient to parameterize the matrix J following Ref [^, using a unitary transformation 
U to bring it into upper triangular form. 



T = UUU 



A2 



T2N 







A 



(33) 



N j 



The ensemble requires 2'N'^ coordinates. Of these, 2A^ are given by real and imaginary parts 
of the eigenvalues Ajt, and N{N — 1) by real and imaginary parts of the matrix elements T^^. 
The remaining N[N—1) parameters Hki are as described by Mehta ||23|. The Jacobian of this 
transformation is proportional to Y[k<i |Afc — A^p and thus depends on Ai, . . . , A^r only. Note 
also that the eigenvector correlator O^/? is invariant under the unitary transformation U. In 
this section, (L^| and \Ru) will denote left and right eigenvectors in the new basis. Thus, 
T\Ri) = Xi\Ri) and \Ri) = (1, 0, . . . , 0)^. In keeping with Eq. (D, let (Li| = (1, 62, • • • , &iv). 
The coefficients bi can be determined by recursion: From {Li\T = AiT one has, with bi = 1 
and for p > I, 



Ai — A, 



2_^ 1 IP 
9=1 



(34) 



The solution of this recursion relation is 
61 = 1 

02 



Ai-A 



2 



, _ ^13 T12T23 , . 

Ai — A3 (Ai — A2)(Ai — A3) 

Trri rri rri rri 

, 14 , -'13-' 34 , -'12-' 24 

^4 = -; ;- + 7^; r^rr^; r^ + 



+ 



A1-A4 (Ai-A3)(Ai-A4) (Ai-A2)(Ai-A4) 

T12T23TS4: 
(Ai-A2)(Ai-A3)(Ai-A4) 



Eq. (pSD provides an explicit expression for the correlator 



N 



On = Y,\bi\' (36) 



1=1 



in terms of the eigenvalues A^ and the matrix elements T^; for k < I. 

To calculate off-diagonal correlators one needs, in addition, the eigenvectors (L2I and |-R2)- 
Let |i?2) = (c, 1, 0, . . . , 0)^ and, in keeping with Eq. (D, (L2I = (0, 1,4,..., dj^). Eq. (|) 
implies that c = —62- Then {L2\T = A2T gives, with di = and d2 = 1, 



d 



V 



1 "'^ 



A2 — Ap 



9=1 



This recursion relation is solved in the same way as (Bl) and 



N 



Oi2 = -h2Y,dih (38) 



1=1 



provides a corresponding expression for O12 in terms of the eigenvalues A; and the matrix 
elements Tki {k <l). 

2. Integration on T^i 

It was shown in the previous section how the correlators On and O12 may be expressed in 
terms of the eigenvalues A^ and the matrix elements Tm {k < I). The Jacobian depends 
only on Ai, . . . , \n- In calculating averages of the the type ([l0|) and (0), the A^(A^ — 1) 
parameters Hki mentioned in section |111B 1| can thus be integrated out and only the integrals 
over Tfci for k < I and over A^ for k = 1, . . . ,N remain. These have the form 



fU'i'x, n 1^' - ^'1' nd'r„ . . . exp f-i y: i^^i' - ^ e \t.A 

•^ k k<l k<l \ k k<l J 



(39) 
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The integrals on all the eigenvalues will be discussed in the next section. In the present 
section we show how to perform the integrals over T^i. To this end the notation (■ ■ ■)j' is 
introduced, denoting a normalized integral on all T^i with weight exp(— o"~^ Sfc<i l^fciP)- 
Consider first the average {Ou)t- Let 



Si = J2(\^p\')t 



p=l 

so that 5*1 = 1 and Sn = (Oii)^- Then from Eq. 

,2 



\bi\'h 



a 



|Ai-A,P 



and hence 

Together with 5*1 = 1 this implies 



a 



-A/P 



Si-i 



Si-i 



N . 
/=2 ^ ' -"- 



a 



A,P 



Consider now the average {0\2)t- Let 

I 
S'l = {T 12^ bkdk)T 



k=l 



SO that S[ = 0,3'^ = aV(Ai - A2) and {Ou)t = -S'j^/{Xi - A2). Now 
Si+i - Si = (T12 bi+i di+i)T 



T 



12 



1 

— Yl blTql + l 



Ai — A/+1 



q=l 



^2 - M+1 fc^i 



T 



(Ai — Ai+i)(A2 — A/+i) 



S[- 



This implies 



(Oi2)^ 



a 



N 



Ai — A2 



n^ 



i=3 



a 



(Ai-A0(A2-A0 



(40) 



(41) 



(42) 



(43) 



(44) 



(45) 



(46) 



Eqs. (1^) and (|46|) represent the averages of On and O12 with respect to the coordinates 
Tfc/. The remaining integrals are those over A^. Using Eqs. (^31) and (E3) one has 



N 

1=2 ' ^ 



a 



\i\y/p 



(47) 
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and 



0(zi, Z2) = -N{N - 1) a^ {6{zi - Ai) 6{z2 - A2) 



a 



N 



X 



Ai — A2 



n 1+ 



a 



where 



,^3 (Ai-A0(A2-Az) 

■)p is an average with the weight (PTI). 



(48) 



3. The case N 



The case A^ = 2 is particularly simple. We find 



0{z) = -[2 



and 



0(^1,2:2) 



^2^2 



a^ 



exp 



exp 



a^ 



Fil 



F2I 



a^ 



a^ 



(49) 



(50) 



These expressions are useful as simple checks of results for arbitrary values of A^. 
Of more general interest, the distribution of Oaa is, from Eqs. (p^) and (p6D 



P{Oc 



QjOaa - 1) 
(2 0, 



1)3' 



where B(x) = 1 for x > and zero otherwise. This gives 

(O,,) = 3/2 



(51) 



(52) 



which is consistent with (^9D integrated over z. Note that the second and higher moments 
of Oaa diverge. We argue in section |V A2| that the the tail of the distribution of Oaa at 
large Oaa has the same form for all A^ > 2. 



4- Calculation of the eigenvalue averages 

In this section we show how to evaluate the remaining integrals in (^) and (^8]). They can 
be performed in the same way as those in section |111 A| . In analogy with Eq. (^) one has 



0(z) =Na^(7TN\a 



2\-l 



exp 



a^ 



det 



900 901 
910 







9N-3N-2 
9N-2N-3 9N-2N-2 . 



(53) 



withgij = {TTj\a^^+^y^Jd^XVx'{a^+\z-X\^) exp(-a-2|An. Eq. (H) provides an explicit 
expression for 0{z) for general A^. The determinant can be easily evaluated numerically, as 



12 



is shown in section 0. For z = 0, the (A^ — 1) x (A^ — 1) determinant in Eq. ( p4D is simply 
diagonal. Denoting it by Gat-i we have G^+i = {k + 2)Gk and thus 



Oiz] 



_ N 
z = n 



(54) 



independent of cr^. For N = 2 this expression gives 0(0) = 2/7r, consistent with Eq. ( ^9D 
An expression for 0{zi, Z2) can be obtained in analogy with (p8|) : 



0(zi,Z2) = -N(r^n'N\a 



,4\-l 



exp 



I 12 I 12 

\Zl\ \Z2\ 



a^ 



a^ 



(55) 



X det 



^00 ^01 ^02 
hiQ hii 

h2Q 



hN-5N-3 
hN~3N-5 h^-ZN-i hN-3N-3 



with 



+a\z^~X){z2-X) 



Z^~X\^\Z2-X\' 



(56) 



exp 



a^ 



For Z2 = and zi = z, and denoting the determinant in (^) by Hn-2i we obtain the 
recursion relation 



Hk+i = {ex-' \z\' + k + 3)Hk - a'' \z\' {k + 2) Hk-i . 
With Hi = 3 + (T-2 l^p and H2 = 12 + Aa-^ \z\'' + a"^ \z\^ this yields 

|2s 



O(z,0) = -Ar.^vr-exp(-^)^5: 



^ 'zP' 



a^ / |zr ■' — ' /! a 
1=2 



21 



(57) 



(58) 



For N = 2 this gives O{z,0) = — (crvr)"^ exp(— cr~^ kP), which is consistent with (|50|). An 
additional check is provided by the fact that Eqs. ( ^4|) and (|58D obey the sum rule (|l6l). In 
the limit of A^ large, with cr^ = N~^, we obtain, for z y^ 0, 



O{z,0) 



7t'\z\^ 



(59) 



for l^l < 1 and zero otherwise. In order to exhibit the behavior of Eq. ( p5| ) near the origin, 
for a^ = N~'^ and in the large A^ limit, we write u = N^^'^z; for \uj\ -C N^'^ we then have 



A^"^O(2,0) 



7r2|cj|4 



l-(l + |a;P)e-H^ 



(60) 



Eq. (pOD displays the way in which the result (|59D is regularized as \z\ — > 
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5. Simplified calculation of the eigenvalue averages for N large 
The main results of section [111 B 4| are the determinantal expressions Eqs. (|5B|) and 



providing exact results for the eigenvector correlators ([T0|) and (|TT]). In the present section 
we provide approximate expressions for ( ^31) and (|55|) which, for a^ = A^~^, are valid in the 
limit of A^ large, with zi ^ z^ and |zi|, |z2| < 1- In the following we shall need to indicate 
explicitly the rank of the random matrix considered, and so we use the notation 

in place of 0(z) and 0(^1,^2) [Eqs. ( p!OD and ([11]) ]. Consider first Ojq{z). We write 



0^(;2i, z^) = l^Y. ^-^ -^(^i - ^") ^(^2 - A^)) . (62) 



Om{z)^Om{^)^i (63) 

where 

N 

".- n (i+]v^3i^) (*'^' 

and the product excludes the M eigenvalues A; closest in the complex plane to the point z, 
as illustrated in Fig 1. We believe that Eq. (|6^) is exact for A^ -^ cxo followed by M -^ 00, 
because we expect that Vti has no fluctuations in that limit. This implies in particular that 
we can calculate Vti by evaluating the average of its logarithm. Starting from 

N 

logn.= 5:iog(l + ^^— ^) (66) 

l=M+l ' '' 

and expanding the logarithm on the right side, we have 

l=M+l ' i.\ -J U I I 

where, in the large A^ limit, d{z) = vr"^ for \z\ < 1 and d{z) = otherwise. The domain D 
of integration excludes a disk of radius g with centred on z. Since this disk should contain 
M eigenvalues, g^ = M/N. Thus we obtain in the large N limit 

Q^ = g-' (1 - \z\') . (67) 

Making use of the fact that Oa/ (0) = M/n [see Eq. (pif)], and using Eq. (|63|) we thus obtain 

0{z) c:, —{1 - \z\^) . (68) 

IT 

The quantity Oi\f{zi, Z2) can be calculated in a similar fashion. To this end we write 
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On{zi, Z2) = Om{zi - Z2, 0)^2 (69) 

where fl2 is 

N 

n2= Y\ (i + ^T ^) (70) 

and the product excludes the M eigenvalues closest to Z2, with 1 -C M <^ A^. We first 
consider the case \zi — Z2\ ^ I ■ Proceeding as above, we have 

log n2= [ d'x d{x) l_ (71) 

where again the domain of integration T> is the unit disk with a disk of radius g around Z2 
removed as illustrated in Fig. |I[ In the large A^ limit we obtain 

^2 ^ f?"'(l - Z1Z2) . (72) 



As before, ^^ = M/N. Using Eqs. (58) and (^), we find in the large A^ limit and with 

uj = Ny^Zi-Z2), \uj\ <Ar^/2 

N-'0{z,, Z2) ^ - ^^ (1 - (1 + \u:\') e-H^) . (73) 

Second we consider the case \zi — ^2! ^ N"^^"^. In this case we obtain, in the large A^ limit, 

for Zi 7^ Z2 and l^il, 1^2 1 < 1, 

Oh,..) = -ii^. (74) 

vr^ \zi — Z2\ 

For \zi\, \z2\ > 1, 0(21, 2:2) vanishes in this limit. 

IV. GIRKO'S ENSEMBLE 

In this section, we present a general approach to calculating the averages of Eqs. (p!o|) and 
(p!T|) perturbatively, using an expansion in powers of A^~^. This will enable us to treat 
more general ensembles than the one considered in the previous section. As an example, 
expressions are derived for the averages of Eqs. (p!0| ) and ( pJ]) in the case of Girko's ensemble. 



defined in Eq. (^) . The expressions derived below are appropriate for large A^ and zi 7^ Z2 in 
(|n|). For r = 0, Eqs. (^^ and (^) are thus reproduced. In the following we set a = N~^^'^: 
the results derived are correct in the large A^ limit. 

A. Self-consistent Born approximation 

The desired approximations for ( PUJ ) and ([TT|) are obtained by calculating the average in 
Eq. (pOD using Green functions. The corresponding Green functions are non-analytic within 
the support of the density of states which occupies a finite region in the complex plane. 
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In general, perturbation theory yields only the analytic contribution, and in conventional 
problems singularities on the real axis are obtained by analytic continuation. In the present 
case one thus proceeds as follows. A Hermitian 2N x 2N matrix H = Hq + Hi is introduced 
BJlJTlJTlilJll 



Hn 



-7] 



H^ 



A 



A^ 



with ?7>0,y4 = z — J and with inverse 



G 



A^[rf + AA^]-^ -'nlrf + A^ A]-^ 



Gil Gi2 
G21 G22 



(75) 



(76) 



Expanding the Green function as a power series in Hi, its ensemble average {G) can be 
written as 



(G)=Go + GoS(G) 



(77) 



where Gq = Hq ^ and S is a self-energy. Within the self-consistent Born approximation one 
obtains ||7|JT5| 



IN 



{G22) ~z + r{G2i) 
-z + t{Gi2) {Gil) 



(78) 



as illustrated diagrammatically in Fig. |^. The self-consistent Born approximation is exact 
in the limit N ^ 00. For t] ^ 0, the self-consistent solution of Eqs. (|77D and ([78|) is as 
follows P,|7|,p!5[] : one has for all z {G22) = — (Gii) and (Gu) = {G2i)- In addition, (Gn) is 
non-zero only inside the ellipse defined by [x/{l + r)]^ + [y/{l — r)]^ = 1 



(Gil) 



1 — [x/{l + r)]2 — [y/{l — r)]2 inside the ellipse 



Furthermore, 



(^21) 



x/(l + r) - iy/ {l - t) 
{z - V^2 _ At) /2t 

Using Eq. ([ISD , the density of states is given by 

1 d 
d{z) = lim-7— {G21) . 
v^o ndz 



outside . 



inside the ellipse 
outside . 



(79) 



(80) 



(81) 



It thus turns out that for A^ :^ 1 the support of d{z) is an ellipse in the complex plane 
iJ^jTl with 

^.,^ ^1^-^(1- r')-' for [x/(l + r)]^ + [y/(l - r)]^ < 1 , 
\0 otherwise . 

In the limit r ^ 1, the eigenvalue density of the Gaussian Unitary Ensemble is recovered, for 
which d{z) = 6{y) (27r)"^-\/4 — x^. Alternatively, setting r = 0, the support of the density of 
states in the complex plane becomes a disk of unit radius centred around the origin [compare 
Eq. 
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B. Bethe-Salpeter equation 

In the following, Gki{zi,zi) is denoted by Gki{l) {k,l = 1,2). An equation for the average 
of the matrix product (G2i(l) ^12(2)), accurate at leading order in N~^, is shown diagram- 
matically in Fig. ||. There are sixteen such equations for all products {Gij{l) Gkii^)) for 
i / = 1, 2. In order to write these in matrix form, one defines 



i?(l,2) = (G(l)®G(2)). 



(83) 



Similarly, Rq{1,2) is the matrix (G(l)) ® (G(2)). The matrices R and i^o are Hermitian. 
Defining the vertex 



/ 



In 



T 



T 



Vi 



1\ 



/ 



(84) 



the diagrammatic expression for R{1, 2) can be written as 

il(l, 2) = Roil, 2) + Roil, 2) r Ril, 2) . 
Eq. (|86D has the solution 

Ril,2) = [l-Roil,2)T]-' Roil,2). 



(85) 



(86) 



We first discuss the simplest case, r = 0. If 2:1 and Z2 lie inside the support of the density 
of states i\zi\ < 1 and \z2\ < 1) 



IN 



i?o(l,2) 
In this case, from Eq. (|86|) 

/ yi - ZiZiJl - ZiZ2 



Z\ 



vT^^iZi 

zi -VI - 2:1^1 



® 



Z2 



VI - Z2Z2 

Z2 -VI - Z2Z2 



(87) 



R(l,2) 



IN'- 



'1 - ZiZ 



'1 - ZiZi 1 - ZiZi - Z2Z2 + Z1Z2 \ 



\Zl - Z2\ 



_ y^l - ZxZx 

yi - Z2'Z2 
Z\ - Z2 

1 — Z\Z\ — Z2Z2 + Z\Z2 
\Z\ - Z2? 



Z\ - Z2 



^1 -^2 



^1 - 2^2 



^1 - Z2 

Z\ — Z2 

^1 -Z2 





\Zl - Z2\ 
yi - Z2'Z2 
^1 -2^2 

/I - 2:1^1 
Zl - Z2 



'I-Z2Z2 Jl-ZiZi Jl- ZiZiJl- Z2Z2 



Zl - Z2 



Zl - Z2 



\Z\ - Z2\ 



(88) 



Alternatively, if both \z\\ > 1 and \z2\ > 1, we obtain 



17 



R(h2) 



IN 



1 \ 

[Z2 - 1 



2:1^2 



Z1Z2 



(89) 



/ 



The general case, r 7^ 0, is dealt with as follows. We define a transformation 

Z — TZ 
Z ^' W 



l-r2 



(90) 



which maps the support of the density of states in the 2;-plane onto the unit disk in the 
w-plane. For zi and Z2 inside the support of the density of states one has \wi\ < 1, \w2\ < 1 
and 



Ro{l,2) = l 



N 



y/l — WiWi Wi 

Wi — VI — WiWi 



VI - W2W2 W2 

W2 -a/1 - W2W2 



The resulting matrix -R(l, 2) is more complicated than Eq. 
(G'2i(l)G'i2(2)) in the case \wi\ < 1 and \w2\ < 1 we find 



(91) 
For the element 



(G2i(l)Gi2(2)) (92) 

_ (1 - r)^ + (1 + t'^)uJiW2 - WiWi - W2W2 + t{wiW2 + W1W2 -w\- wl) 

IWi - W2 + t(Wi - W2) P 



C. Calculation of the density D{zi,Z2) 

The density D{zi, Z2) can be expressed in terms of Green functions, from Eq. 

D{z,,Z2) = hm- — — (G2i(l)Gi2(2)) . 

»?^0 TT^ OZi OZ2 

We find from Eq. (p^), for \zi\, \z2\ within the ellipse, that 

(1 _ rf (1 - T^f - (1 + T^)z{Z2 + t{zI + ^l) 



D{Zi,Z2) 



TT^ 



Fl - Z2 



as 



(93) 



(94) 



For zi and Z2 outside the ellipse, D{zi, Z2) vanishes. 

As a check it can be shown explicitly that D{zi, Z2) obeys the sum rule (p!6|). Using Green's 

theorem, we have 



fd'z2D{z,,Z2) = ^ J^^ /dZ2(G2l(l)Gi2(2)) 



(95) 



where the contour integral is around the ellipse. By means of the transformation (pOD, this 
contour may be mapped into the unit circle in the w-plane, giving 



|d^-,flfe,.,) = iT^(^-r^)i/(d-. + r*.)(G.,(l)G..(2)> 



7rl-r2 



|«'2| = l 

(96) 



for \wi\ < 1 and zero otherwise, as expected from Eq. (82). 



As a final check we observe that, with r = 0, Eq. (M) imphes 



0(^1,^2) = -^^^ (97) 

TT^ \Zi - 2:2 f* 



for \zi\, 1^2! < 1 and zero otherwise. Thus, our previous result, Eq. (ff^) , is reproduced from 
(0) for r = 0. 

As pointed out in section |TB|, the diagonal correlator 0{z) is given in terms of the singular 
part of D{zi,Z2), see Eq. (P^l). This singular part is inaccessible perturbatively, in lowest 
order in A^~^ [^. In order to determine 0{z) within the perturbative approach discussed 
in this section, we proceed as follows. For simplicity, consider the case r = 0. Integrating 
the density D{zi, Z2) over a small disk around Z2, of radius r] which is taken to be small 

j d'z2 D{z,, Z2) = ^^j dz2 ^ J- (G2i(l) Gi2(2)) ~ ^(1 - IziH , (98) 

\zi-Z2\<r) zi-Z2=ri 

provided Zi is sufficiently far away from the boundary. On the other hand, from Eq. (pjlf) 
and for r] ^ A^^^/^, this is approximately 0{zi), so that up to prefactors of order 0(1), 

0{zi)^N{l-\zi\'^) (99) 

[compare Eq. (|68|) ] and thus Oaa ~ A^- The sum rule (|^) can be used to check the 
consistency of Eqs. (0) and (|99D. 



V. SUMMARY AND DISCUSSION OF THE RESULTS 

In the present section we summarize and discuss the results obtained in the previous two 
sections. As in Sec. |^, the variance a^ in Eqs. (|I]) and (0) is taken to be 1/A^. 

A. Ginibre's ensemble 



1. Eigenvector correlators Eqs. ^lUj ) and JT^ j 



In the case of Ginibre's ensemble we have been able to obtain exact expressions for the 
eigenvector correlators, Eqs. (0) and (|ll]), in the form of determinants. In certain cases, 
we could simplify these expressions further by recursion. Combining these results [compare 
Eqs. (|5^) and (^8])] with a continuum treatment (see section |IIIB5| ), in a way which we 
believe gives exact results for the large A^ limit, we have for |-2i — -22! 7^ and \zi\, \z2\ < 1 
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Ar-io(^i) = -(l-|^i|2) (100) 

0(^1,^2) = -^^^^ (101) 

TT^ \Zi - Z2Y 

For \z\\^ \z2\ > 1, both densities vanish as A^ — >■ 00. To display the form of 0{zi,Z2) as 
\zi — Z2\ — > 0, it is necessary to express zi — Z2 in units of the separation between adjacent 
eigenvalues. Let z+ = {zi + Z2)/2, Z- = z\ — Z2-, and u) = \fNz-. For \zj^\ < 1, a; <^ \/N 
and for A^ ^ 1, Eq. (|73D implies 

N-'0{z,, Z2) = -^^ (1 - (1 + W) e-l^l') . (102) 

We have examined the convergence towards these results for increasing N . In Fig. ^ we show 
A^~^ 0{z) as a function of z for A^ = 2, 4, 8 and 16, obtained by evaluating the determinant in 
(0). We also compare this with Eq. (|100|) . The exact results converge rapidly towards the 



approximate result ( [lOOj ) as A^ is increased, provided z is sufficiently far from the boundary 
of the support of 0{z). 

In Fig. |] we show 0(zi, Z2) as a function of Zi (on the real axis) for Z2 = 0.4 for A^ = 2, 4, 8 
and 16, obtained by evaluating the determinant in (|55|). We compare this with Eq. ( |101[ ). 
Again, the exact results converge rapidly towards the approximate expression ( |101| ) as A^ 
is increased, provided zi and Z2 are not too close to each other or to the boundary of 
the support of 0(21,^2)- Finally, in Fig. |] we show the behavior of N~^0{zi,Z2) for 
\zi — 2:2! ~ A^~^/^, comparing the approximate expression (|102D with exact results obtained 



by evaluating the determinant in Eq. (|55D. The exact results converge very rapidly to the 
approximate expression as A^ is increased, provided z+ < 1 and |a;| <C vN. 
It is important to stress the dramatic difference between the behaviour of Oai3 in Ginibre's 
ensemble and its behaviour in the case of Hermitian matrices, for which Oq,^ = 6af3- The 
fact that, by contrast, Oaa ~ A^ in the non-Hermitian ensemble can be understood as the 
behaviour which would result if {L^l and \Ra) were independent random vectors, subject 
to the normalisation of Eq. (|^): Choosing a basis and scaling in which \Ra) = (1, 0, ... , 0)"^, 
and assuming that (Lq| is a random vector, biorthogonality requires {La\ = (1, &2, • • • , ^at), 
where the coefficients bj, for j > 1 are random and \bj\ is expected to be of order 0(1). Thus 
(Oaa) ~ N. Moreover, large values for the diagonal elements of the matrix 0^/3 must be 
accompanied by some large (or many small) off-diagonal elements, since the two are linked 
by the sum rule (||). Indeed, Eq. (0) implies 

Oo,pr^O{z^,Z2)/R2{z^,Z2) (103) 

and hence, from Eq. ( p.02| ), Oap ~ — A^ if Aq and Xp are neighbouring eigenvalues in the 
complex plane, so that (typically) a; ~ 1. 



2. Distributions of O, 



a/3 



Finally, it is interesting to ask about, not only the average behaviour of the overlap matrix, 
but also its fluctuations. In fact, Oap is typically large if the matrix J has an eigenvalue 
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which is almost degenerate with A^ or A/3, and as a result, the probability distribution of 
Oq/3 has a power-law tail extending to large \0ai3\- To illustrate this, we consider N = 2, 
for which the probability distribution, P{Oaa), of a diagonal element of the overlap matrix 
is given by Eq. (|5l|) and decays at large Oaa according to P{Oaa) ~ 0~^. This implies in 
particular that the second and higher moments of Oaa diverge. 

For N > 2, the tail of the distribution P{Oaa) is determined by pairs of eigenvectors with 
closest eigenvalues, and we expect that for general A^, the tail of the distribution function 
decays algebraically according to 

P{Oaa) ~ O-'^ . (104) 

In Fig. 1^ we show the distribution P{Oaa) of the diagonal overlaps Oaa in Ginibre's ensemble 
for A^ = 10. The tail of the distribution function is well described by Eq. ( |104[ ). 



B. Girko's ensemble 



The main result of section |V|is Eq. (0), giving 0(2:1, Z2) provided \zi — Z2\ is much greater 
than the mean separation in the complex plane between neighboring eigenvalues. For r = 0, 
Girko's ensemble reduces to Ginibre's ensemble. Correspondingly, the perturbative result 



reproduces, for r = 0, \zi — Z2\ 7^ 0, \zi\, \z2\ < 1 and large A^ the expression (|101D 



which was obtained from the exact results of section |T| in the same limits. The singular 
contribution of the diagonal overlap matrix elements to D{zi, Z2) is only indirectly available 
within perturbation theory. Eq. (|99|) shows that the singular behaviour extracted from 
the perturbative results is consistent with the exact expressions [compare Eq. (|100|) ]. On 



the other hand, Eq. ( plj) implies that for 1 — r ^ 1, Oa/B ~ 0{zi, Z2) / R2{zi, Z2) oc 1 — r. 
Thus Oai3 vanishes in the Hermitian limit r — > 1, as expected. The same is true for the 
ant i- Hermit ian limit, r ^ — 1. 



VI. IMPLICATIONS 

Fluctuations of eigenvectors in non-Hermitian random matrix ensembles exhibit a num- 
ber of striking features which are likely to be relevant in physical applications. As in the 
immediately preceding sections, we take the variance in Eqs. (|l|) and (0) to be a^ = A^^^. 

A. Sensitivity to perturbations 

First, as pointed out in the introduction, systems described by a non-Hermitian operator are 
particularly sensitive to perturbations. This sensitivity is determined by the diagonal matrix 
elements of Oai3- In order to illustrate this fact, it is convenient to consider a one-parameter 
family of matrices 

J= Ji cose + J2 sin^, (105) 

where the parameter 6 is real and the matrices Ji and J2 are drawn independently from the 
same ensemble. Then 
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{\dXo^/de\^)=N-^0^^). (106) 

According to Eq. ( |10CI| ), {Oaa) is large, being of order N. Thus {\d\a/d6\'^) is of order unity. 
This should be compared with the Hermitian case E^, where {\d\a/d9\^) is of order N~'^ 



in the corresponding parametrization. Structural stability, on the other hand, requires that 
the level velocities tend to zero as the boundary of the support of the density of states is 
approached: the latter must remain unchanged as 9 varies, since the perturbations merely 
take J from one realisation of the ensemble to another. The expression (|100|) for 0{z) shows 
that this is indeed the case. 



B. Time evolution 

Systems governed by a non-Hermitian evolution operator may exhibit transient features in 
the time-dependence of correlation functions which are controlled by the type of correlations 
between left and right eigenvectors that we have studied. Consider for example an evolution 
equation of the form 

^K) = (J-l)|w,) (107) 



with J drawn from Ginibre's ensemble. We use J — 1 rather than J in Eq. ( |107| ) for 
convenience, to suppress exponential growth. This corresponds to shifting the support of 
the density of states by unity along the negative real axis, so that all (except a vanishing 
fraction) of the eigenvalues have negative real parts. Then 

\ut) = Y,\Ra) ft{K) {L^\uo) , (108) 

a 

with /j(A) = exp[(A — l)t]. Ensemble averaging with {uq\uq) = 1 yields 



{{nt\ut)) = (^ E^"/3 /*(^") fti^P)) ■ (109) 



af3 



Thus, properties of the matrix Oap directly influence time evolution. Eq. ( |109|) can be 



obtained as the double Laplace transform of the density (0), with respect to zi and Z2, 

{{ut\ut)) = fd'zid^Z2e^'^+^'-^^'D{zi,Z2) . (110) 

The diagonal and non-diagonal contributions to D{zi,Z2) yield large contributions to Eq. 
( p.lO| ) which almost cancel. It is thus convenient to evaluate the double Laplace transform 
110| ) by contour integration. 



{{ut\ut)) = -^ Lz^Lz2 e(^^+^^-2)* (G2i(l)Gi2(2)) . (Ill) 



\zi\=l \Z2\ = 1 

In this case one obtains for large N and for t <C vN 
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{{ut\ut)) = e-^' Io{2t) (112) 

which, for 1 -C t -C vN, simphfies to 

(113) 



VAnt 

This behaviour should be compared with the much faster decay that would result from the 
same spectrum if the eigenvectors were orthogonal. In the same regime, the replacement 
Oaf3 —>■ Sai3 transforms Eq. (|109D into 

a 

and thus 

{{ut\ut)) = t-'e-''h{2t) ~ ^^ip- (115) 

Thus, eigenvector correlations may be as significant as eigenvalue distributions in deter- 
mining evolution at intermediate times, a fact of established importance in hydrodynamic 
stability theory [|l|,|g]. 



The more general case of Girko's ensemble (|^) can also be treated in this way. Mapping the 
corresponding contour integrals to the u7-plane by means of (pO|), we obtain 

{{u,\ut)) = e-2(i+-)* |/o[2(l + T)t] +Th[2{l + r)t]]. (116) 

Here the support of the spectrum was shifted by 1 + r along the negative real axis. The 
limiting cases of Girko's ensemble, r — >■ ±1, are easily understood: In the anti-Hermitian 
case, for r = —1, one has simply {{ut\ut)) = 1 because all eigenvalues have vanishing 
real parts. In the Gaussian unitary ensemble, for r = 1, on the other hand, for large N, 
d{E) = (27r)-V4-^2 for \E\ < 2 and thus 

{{ut\ut)) = ^ [ t/E v/4^rB^e2(^-2)* = {2t)-'e-^'h{4t) , (117) 

which corresponds to ( p.l6|) for r = 1. For the three cases, r = —1,0 and 1, {{ut\ut)) is 
shown in Fig. ^ (full lines) together with the corresponding asymptotic expressions (dashed 
lines) valid for t :^ 1. 

C. Correlations of eigenvector components 

The space-dependence of correlation functions of more general ensembles (such as the one 
discussed in Refs. ^ and |P) can be modelled by correlation functions of the components 
of {La\ and l-R/?). Under a change of basis given by a unitary matrix U, the components of 
say \Rj3) transform according to (Jl-R/3) -^ {j\U\Rf3) = J2m^j'm{''^\Rf3)- Correspondingly, 
{La\i) — > {La\U^i) = J2i{La\l)Uii. Duc to the invariance of the ensemble under unitary 
transformations, we can write 
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{{La\t){j\Rf,)) = J2{UjmUu)u{La\l){m\Rp) . (118) 



ml 



where (■ ■ ■)u denotes an average over the unitary matrices U. With {UjmUii)u = N~^6ij6im 
this imphes immediately {{La\i){i\Rp)) = N^^Sa/sSij. Consider now averages involving 
four eigenvector components. The only non-vanishing (and non-trivial) averages which are 
invariant under the scale transformation (|^ are 

{{^\R^){L^\J){J\L^){R^\^)) = ^^^(% + {O.f,)) - ^]^^(1 + ^^AOafs)) (119) 



and 



{{t\R^){L^\t){j\Lf,){Rf,\j)) = ^^^(1 + 6..A0a^)) - ^]vrn('^^.- + (Oaf^)) ■ (120) 



Summing Eqs. ( |119D and ( |12(J| ) over i and j one obtains (Oajs) and unity, respectively, as 



expected. It should be noted that the dependence on i,j and a,P does not necessarily 
factorize. This is likely to be of importance in problems with spatial structure. The above 
considerations show that interesting space dependence of correlation functions may arise 
from non-Hermiticity. 

VII. CONCLUSIONS 

In this paper we have analyzed correlations of eigenvectors in non-Hermitian random matrix 
ensembles. Such correlations are of interest partly because they determine some aspects of 
the behavior of systems represented by non-Hermitian operators: for example, such systems 
are particularly sensitive to external perturbations and correlation functions may exhibit 
transient features in their time dependences. As emphasized in the introduction, there are 
numerous instances in which random non-Hermitian operators appear in the description 
of physical problems, and we hope that the results and methods summarized here will be 
of interest in a number of contexts. In particular, we have obtained the following results. 
We have characterized exactly the eigenvector correlations in Ginibre's ensemble of non- 
Hermitian random matrices. We have shown that the sensitivity of the eigenvalues with 
respect to external perturbations is larger by a factor of A^ (where A^ is the rank of the matrix) 
than the equivalent for Dyson's ensembles of Hermitian matrices. Moreover, we have shown 
that eigenvectors associated with two different eigenvalues exhibit strong correlations which 
decrease algebraically with increasing separation between the eigenvalues in the complex 
plane. We have also shown that the probability distribution function of eigenvector overlaps 
has algebraic tails. This implies that fluctuations are large, in the sense that higher moments 
of the eigenvector overlaps diverge. In addition to exact calculations specific to Ginibre's 
ensemble, an alternative, perturbative approach has been developed and used to derive 
corresponding results, in an approximate way, in Girko's more general ensemble of non- 
Hermitian random matrices. In the appropriate cases and in the limit of large A^, the exact 
results are reproduced. 
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FIGURES 




FIG. 1. Eigenvalue distribution in the complex plane, as discussed in Sec. [II B 5 



27 



(a) Diagrammatic notation for the ensemble (0): 



(Gii(2;)) =-K>^, {G2i{z)) = ^^^^ , 



(^12(2;)) =-M|^, {G22{z)) = ^®^ ■ 



\JklJkl) = -K*= =»+- , \JklJlk) = -►*= 

k I Ik k I 



(b) Self energy S: 



/ k 



-*^M^^^- 



-M^^=>o;-+-i- -*->-^>^ 



FIG. 2. (a) Diagrammatic notation for the ensemble (y). (b) Self-energy S in equation (|77|), 
to lowest order in N~^ ^ compare Eq. ([7^). Note that a summation over internal indices in closed 
loops incurs an additional factor of N . 






*^>^" 
^m^. 



*^>«*{ 



=^#**^' 



-t^*^ 



FIG. 3. Equation for the matrix product (G2i(l) Gi2(2)). The diagrammatic rules are as in 
Fig. |. 
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0.5 




FIG. 4. N~^ 0{z) = (A^^^ ^^ Oqo (5(2; — Ai)) as a function of z = rexpit/?. The ensemble 
average is independent of ip. Results for A^ = 2, 4, 8, 16 are shown, together with Eq. ( |100| ) valid 
for large N ( ). 




FIG. 5. 0{zi, Z2) = {N ^ J2a^i3 Oap 5{z2 — Xa) (5(^2 — Xp)) for Z2 = 0.4 as a function of zi on 
the real axis. Results for N = 2,4, 8, 16 are shown, together with Eq. ( |101j ) valid for large N 

( )• 
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N=2 



N=4 
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-0.04 - 
-0.06 





-0.06 
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N=16 





N=32 



N=64 





FIG. 6. Shows Eq. ( |102| ) for z+ = 0.4 as a function of uj for several values of N (- 
shown are the corresponding exact results ( ). 



). Also 
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FIG. 7. A histogram of P{Oaa) of the diagonal overlaps as a function of Oaa/^ in Ginibre's 
ensemble for A^ = 10 ( ). Also shown is the theoretical estimate for the tail of the distribution 

( )• 
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FIG. 8. Shows {{ut\ut)) as a function of t for Girko's ensemble ( ), for r = 1,0 and —1. 

Shown are Eqs. ( |112| ) for r = and ( [LITD for r = 1. For r = —1, one has just {{ut\ut)) = 1. Also 
shown are the asymptotic expressions valid for t ^ 1 ( ), compare Eq. (^Ac 



31 



